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ABSTRACT 

According to popular progenitor models of gamma-ray bursts, twin jets should be launched by the central engine, 
with a forward jet moving toward the observer and a receding jet (or the counter jet) moving backwardly. However, 
in calculating the afterglows, usually only the emission from the forward jet is considered. Here we present a detailed 
numerical study on the afterglow from the receding jet. Our calculation is based on a generic dynamical description, 
and includes some delicate ingredients such as the effect of the equal arrival time surface. It is found that the emission 
from the receding jet is generally rather weak. In radio bands, it usually peaks at a time of t > 1000 d, with the peak 
flux nearly 4 orders of magnitude lower than the peak flux of the forward jet. Also, it usually manifests as a short 
plateau in the total afterglow light curve, but not as an obvious rebrightening as once expected. In optical bands, the 
contribution from the receding jet is even weaker, with the peak flux being ~ 23 magnitudes lower than the peak flux 
of the forward jet. We thus argue that the emission from the receding jet is very difficult to detect. However, in some 
special cases, i.e., when the circum-burst medium density is very high, or if the parameters of the receding jet is quite 
different from those of the forward jet, the emission from the receding jet can be significantly enhanced and may still 
emerge as a marked rebrightening. We suggest that the search for receding jet emission should mostly concentrate on 
nearby gamma-ray bursts, and the observation campaign should last for at least several hundred days for each event. 

Key words, gamma rays: bursts — ISM: jets and outfiows — stars: neutron 
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1. Introduction 

Thanks to the discovery of X-ray, optical and radio af- 
terglows of gamma-ray bursts (GRBs), it is now clear that 
most GRBs are situated at cosmological distances (Costa 
et al. 1997; van Paradijs et al. 1997; Frail et al. 1997). A lot 
of progresses have been achieved during the past decade 
(Piran 2004; Meszaros 2006). Especially, through the de- 
tection of GRB 030329, the association of long GRBs with 
supernovae is firmly established (Hjorth et al. 2003), which 
strongly supports the coUapsar model as the energy mecha- 
nism for long GRBs (Woosley 1993; MacFadyen & Woosley 
1999). Theoretically, the collapse of a massive star will most 
likely give birth to a black hole, surrounded by a temporal 
accretion disk. It is a common sense that the accretion sys- 
tem will produce double-sided jets (MacFadyen & Woosley 
1999; Aloy et al. 2000; Rhoads 1999; Meszaros 2002). The 
GRB can be observed only when our line of sight is right 
on one of the two jets. The collimation of GRB ejecta can 
be tested observationally, through various beaming effects, 
such as the achromatic break in GRB afterglow light curves 
(Sari et al. 1999; Liang et al. 2008), the polarization in 
both the main burst phase and the afterglow phase (Lazzati 
2006), the predicted existence of orphan afterglows (Rhoads 
1997; Huang et al. 2002; Granot & Loeb 2003), and the en- 
ergy crisis already noted in some GRBs (Frail et al. 2001). 
In fact, more and more observational evidences have been 
accumulated today, supporting the idea that many GRB 
ejecta might be highly coUimated. 

Current studies on the beaming effects are mostly con- 
centrated on the emission from the forward jet, i.e., the jet 
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moving toward the observer. The emission from the reced- 
ing jet (or the counter jet) is generally omitted. It is inter- 
esting to note that this ingredient recently has been studied 
by a few authors (Granot & Loeb 2003; Li & Song 2004). By 
some simple analytical derivations, Li & Song (2004) argued 
that the emission from the receding jet can be detected in 
a few cases in the non-relativistic phase of GRB afterglows. 
However, previous studies did not consider some important 
effects, such as the action of the equal arrival time surface 
(EATS). Recently, Zhang & MacFadyen (2009) presented a 
two-dimensional simulation of GRB outfiow. The emission 
from the receding jet has also been included in their calcu- 
lations, but they did not investigate the effects of various 
parameters on the receding jet component. 

In this paper, we will present our detailed numerical in- 
vestigation on the emission from the receding jet of GRB 
in the deep Newtonian stage. Although the GRB jet may 
be complicatedly structured (Meszaros et al. 1998; Kumar 
& Granot 2003; Huang et al. 2004), and the circum-burst 
environment may be wind medium and even associated 
with some complex density variations (Meszaros et al. 1998; 
Chevalier & Li 2000; Gou et al. 2001; Wu et al. 2004), here 
we will only consider the simplest situation, i.e, the homo- 
geneous double-sided jets expanding into a homogeneous 
interstellar medium, which is favored by some recent fits 
(Huang et al. 2000a; Yost et al. 2003). 

The structure of our paper is organized as follows. §2 
is mainly a review of the dynamics and radiation model 
we used in our calculations. In §3 we present the numerical 
results, together with our tentative explanations. §4 is our 
conclusion and discussion. 
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2. Model Description 

In the afterglow phase, the GRB ejecta expands into 
the interstellar medium (ISM) and is decelerated continu- 
ously, giving rise to a strong external shock. The swept-up 
electrons arc accelerated by the blastwavc, producing the 
afterglow mainly through synchrotron radiation. In radio 
bands, the shell is no longer optically thin, so that the syn- 
chrotron self-absorption should be considered. In our study, 
we will use the simplified dynamical equations suggested by 
Huang et al. (1999, 2000b), which is consistent with the self- 
similar solution of Blandford & McKee (1976) in the ultra- 
relativistic phase, and is consistent with the Sedov solution 
(Sedov 1969) in the non-relativistic phase. The beaming 
effects (Rhoads 1997, 1999) can also be conveniently simu- 
lated in this way. Here, for completeness, we first describe 
the dynamics and the radiation process briefly. 

2.1. Hydrodynamical Evolution 

In our description, t is the photon arrival time mea- 
sured in the lab frame; R is the radial coordinate measured 
in the burst frame relative to the initiation point; m is the 
rest mass of the swept-up medium; 6 is the half-opening 
angle of the ejecta; 7 is the bulk Lorentz factor of the mov- 
ing material; p is the electron distribution index which is 
typically between 2 and 3; n is the number density of ISM; 

and are the energy equipartition factors for electrons 
and the comoving magnetic field. Wc further denote the ini- 
tial values of the rest mass, the isotropic equivalent energy, 
the Lorentz factor and the half-opening angle of the ejecta 
as Mcj,i?o,iso,7o,fi'j, respectively. 

The overall dynamical evolution of the GRB ejecta can 
be depicted by 
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where /? = -^/l — 1/7^, c is the speed of light, and Cg is the 
comoving sound speed, which can be calculated by = 

7(7 - 1)(7 - l)cV [1+7(7 - 1)] with 7 ~ (47 + l)/(37) 
being a reasonable approximation for the adiabatic index. 
In Equation (4), e is the radiative efficiency. In the extreme 
case, e = means adiabatic condition and e = 1 refers 
to highly radiative situation. Note that in realistic case, e 
should evolve gradually from 1 to 0, in about several hours. 

Equations (1) — (4) is a convenient description of GRB 
afterglow dynamics that is applicable in both the initial 
ultra-relativistic phase and the late Newtonian phase. 

2.2. Radiation Process 

Basically, wc assume that the shock-accelerated elec- 
trons follow a power-law distribution according to their en- 
ergies, d7Vg/d7e oc 7~P, However, to ensure that the cal- 
culation in the deep Newtonian phase is correct, we need 
to modify the basic distribution function as dN^/dje oc 



(7e — 1) ^ (Huang & Cheng 2003). The minimum and max- 
imum Lorentz factors of electrons can be calculated as 

7e min = Ce(7 " l)Wp(p - 2)/[TOc(p - 1)] -|- 1 and 7e,max = 

A/67re/ (o-tS') w 10^{B' /1G)-^^^ , where B' is the comov- 
ing magnetic field strength, rUp and me are masses of pro- 
ton and electron, respectively. As usual, we assume that 
the energy ratio of magnetic field with respect to internal 
energy is , so that the energy density of magnetic field is 

B'^iSn) = (7 - 1)"' (77 + 1)(7 - l)nmpC^. 

The cooling of electrons due to synchrotron radiation 
will lead to a steep distribution function above a critical 
Lorentz factor, 7c. The expression for 7c can be derived 
as 7c = Qirmec/ {a'T^B''^t), where ctt is the Thompson 
scattering cross section (Sari et al. 1998). Considering all 
the above ingredients, we finally use the following electron 
distribution function in our calculations (Huang & Cheng 
2003): 
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In the comoving frame, the synchrotron radiation power 
at v' is (Rybicki & Lightman 1979) 
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with F{x) = a;/^ if5/3(fc) dfc being the Bessel func- 
tion and I'c ~ i^^eB' / (47rmcc) being the characteris- 
tic emission frequency (Shu 1991; Longair 1992). To cal- 
culate the radio afterglows, we must consider the syn- 
chrotron self-absorption. The optical depth of synchrotron 
self-absorption can be obtained as 



Vle^B' 



iinin(7.,n,in , 7c) V d7e 7 7e V ^^e / 



where dng/d7e denotes the column density distribution of 
electrons measured in the comoving frame on the line of 
sight; q is the electron power-law distribution index which 
varies from 2 to p + 1 for fast-cooling and from p to p + 1 
for slow-cooling. The synchrotron self-absorption will affect 
the radiation by a reduction-factor (Waxman et al. 1998) 



f{r) 



1 
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Let us define the Doppler-factor as D = [7(1 — 
(Meszaros 2006), where /i = cos© and 6 is the angle be- 
tween the velocity of the emitting material and the line of 
sight. Also we denote the viewing angle as ^obs- Then the 
observed frequency is u = Dv' /{I + z), and the observed 
flux density from a point-like source is 
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where dh is the luminosity distance. FinaUy, we can inte- 
grate the flux density over the EATS (Waxman 1997; Sari 
1998) determined by 



t, 



obs 



dR 



= const. 



(12) 



3. Numerical Results 



In this section, we present our numerical results con- 
cerning the emission from the receding jet. First of all, for 
simplicity, we assume that the twin jets have the same char- 
acteristics, i.e., the same initial energy, opening angle, ini- 
tial Lorcntz factor, and the circum-burst ISM density. We 
also assume that the microphysics shock parameters {p, ^e, 
^1) are the same for the receding and forward blastwaves. 
For convenience, we define a set of parameter values as 
the "standard" condition: n — l/cm"^, i?o,iso = lO^^ergs, 
= 0.1, e = 0, = 0.1, ^2 = 0.01, p = 2.5, 0obs = 0, and 
7o = 300. These values are typical in the study of GRB af- 
terglows. For redshift, we adopt the value of 2; = 0.1 (which 
corresponds to cJl = 454 Mpc according to the popular 
cosmology model, Wright 2006). 

Firstly we illustrate the evolution of the Lorentz factors 
of the twin jets in Fig. 1. Note that the X-axis is observers' 
time. For the observer, the dynamical evolution of the re- 
ceding jet is quite different from that of the forward jet, 
especially in the relativistic phase. We see that in a rather 
long time (f ^ 50 d), 7 of the receding jet remains almost 
constant. This is due to the time delay induced by the long 
distance between the twin jets. It also implies that the emis- 
sion from the receding jet will be very weak in this period, 
since it is highly beamed backwardly. At the observers' time 
of t ~ 340 d, the Lorentz factor of the receding jet is still 
more than 10, while the forward jet's Lorentz factor has 
already decreased to less than 1.1 . 

In Fig. 2, we show some examples of the equal arrival 
time surfaces (EATSes) at three moments. As expected, at 
any particular moment, the typical radius of the surface is 
much larger for the forward jet branch as compared with 
that for the receding jet branch. Also, we notice that the 
curvature of the two branches is quite different. Generally, 
the EATS is much flatter on the receding jet. Another inter- 
esting feature is that the area of the EATS on the forward 
branch is much larger than that of the corresponding re- 
ceding branch. 

Fig. 3 shows the radio and optical afterglow light curves 
under the "standard" condition (thick lines). Here, the 
thick dotted line corresponds to emission from the forward 
jet, the thick dashed line corresponds to emission from the 
receding jet, and the thick sold line is the total light curve. 
Under the "standard" condition, for the forward jet, the 
afterglow light curve (the dotted line) becomes slightly flat- 
tened in the non-relativistic phase. It is consistent with pre- 
vious results in the deep Newtonian phase (Huang & Cheng 
2003). Also it can be seen that the receding jet really can 
contribute a significant portion in the total emission at very 
late stage. The role played by the receding jet is reason- 
ably more important in the radio band than in the optical 
band. However, the dashed component is generally not very 
strong, so that it can only lead to a plateau in the total light 
curve, but not an obvious rebrightening or a marked peak 
as expected by Li & Song (2004). Interestingly, our result 
is consistent with the simulation of Zhang & Macfadyen 



(2009). We believe that the discrepancy between our nu- 
merical result and Li & Song's analytical result mainly 
comes from the effect of the EATS. Below, we will give some 
detailed analyses on this point. Additionally, it should be 
noted that in the radio band, the peak flux of the reced- 
ing component is about 4 orders of magnitude weaker than 
that of the forward component. It essentially means that 
the receding component is very weak, and is very difficult 
to detect. In the optical band, the condition is even more 
awkward. The peak flux of the receding component is about 
23 magnitudes dimmer than that of the forward component 
in R band. Even comparing with the flux of the forward jet 
at the jet break time, it is still 16 — 17 magnitudes weaker. 
So, in optical band, it is even much more difficult to observe 
the receding jet component. 

According to Li & Song (2004), the time when the re- 
ceding jet becomes notably visible (t^ri) is relevant to the 
time when the forward jet enters the non-relativistic phase 
(iNR), i.e., 



^RJ 
f-NR 



^NR 



2rNR 



(13) 



where tnr is the radius of the forward jet at tNR. 
In the standard frame work (Blandford & Mckee 1976; 

Rhoads 1999), the sphere- like phase of a highly coUimatcd 
GRB ejecta ends at the so called jet break time deter- 
mined by 7j = 1/Oj, with the shock radius being rj = 

1/3 

(3i?o.isof^jV [47rnTOpC^] ) . After the sphere-like phase, the 
jet spreads laterally at the co-moving sound speed Cg so 



that we have 7 oc f and tnr 
finally we obtain 



(Rhoads 1999). Then 



^NR 



1 / 3ii^o,i! 



1/3 



2c \ 4nnmpC^ 



.RJ 
^NR 



S^NR. 



(14) 



Adopting the standard values of our parameters. 
Equation (14) yields tNR « 104 d and « 520 d. After 
correcting for the cosmological time dilation {z = 0.1), we 
get the corresponding observers' time of ti = (1 +z) « 
114 d and = (1 + z) t^^^ « 572 d. In fact, in Fig. 2, 
the EATSes for these two moments have been displayed. 
So, according to Li & Song's suggestion, the contribution 
from the receding jet should peak at ^3 « 572 d. In our 
Fig. 3, for the "standard" condition, the peak is postponed 
to tpeak ~ 1140 d for 8.46 GHz, and to tpeak ~ 1700 d in 
R band. So, the EATS effect and the deceleration of the 
external shock can lead to some subtle difference between 
the analytical results and the numerical results. Actually, 
Zhang & MacFadycn's numerical results have clearly shown 
that the observers' time does not equal to the burst frame 
time at tNR (Zhang & MacFadyen 2009). Unfortunately, 
in previous analysises it is usually assumed that these two 
times are equal. 

Another reason that suppresses the rebrightening of the 
receding jet is as follows. According to Li & Song's analy- 
sis, at the observers' time ta, the receding jet should be at 
the radius of tnr. However, from our Fig. 2, we see that 
the typical radius of the EATS at ^3 on the receding jet is 
much smaller than the radius of the forward jet at ti. The 
reason is again due to the EATS effect. It means that the 
receding jet still does not decelerate enough at ^3 (actually, 
the bulk Lorcntz factor is still 3.95), and its emission is 
still mainly directed forwardly (not backwardly toward the 
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observer). Additionally, Fig. 2 shows clearly that the area 
of the receding jet at (corresponding to t^p^ ) is much 
smaller than that of the forward jet at ti (corresponding 
to tNR,)- So, the number of electrons involved in the ra- 
diation process is typically much smaller on the receding 
jet at t^R' ^ compared to that on the forward jet at ^nr- 
Due to the above reasons, the contribution from the reced- 
ing jet is naturally much weaker than that deduced from 
L^'-^t) w L„{t - 4tNR), {t > t^^) (Equation (7) in Li & 
Song (2004) ). 

However, although the receding jet emission is gener- 
ally very weak in our "standard" condition, we guess that 
in some special cases it still can be enhanced. Obviously, 
a denser environment will help to decelerate the jet more 
quickly, thus lead to a smaller fpeak and a higher intensity. 
In Fig. 3, we have also plotted in thin lines our numeri- 
cal results for a double-sided jet that locates in a dense 
circum- burst medium (n = lOOO/cm"^). Note that other pa- 
rameters involved here are the same as the "standard" case. 
Encouragingly, in Fig. 3(a) we see that the peak time of the 
receding jet can be as early as tpeak ~ 150 d, with the peak 
flux as large as a few mJy in radio band (i.e., only sev- 
eral times less than the peak level of the forward jet). In 
Fig. 3(b), the optical contribution from the receding jet is 
still very weak, with the peak flux being about 28"*. 

In Fig. 4, we plot the afterglow light curves in more ra- 
dio and optical/infrared bands. Generally speaking, ipcak is 
about 1140 d in radio bands and is about 1700 d in optical 
bands. Such a difference in the peak time is insignificant, 
considering that the frequency difference between radio and 
optical wavelengths is really huge. We notice that tpoak al- 
most remains the same from radio to X-ray bands in Fig. 7 
of Zhang & MacFadyen (2009). Thus our results are roughly 
consistent with Zhang & MacFadyen's. Another interesting 
conclusion that can be drawn from our Figs. 3 and 4 is 
that at lower frequency, the relative intensity of the reced- 
ing jet component (its peak flux), as compared with the 
peak of the forward jet component, becomes stronger. Such 
a tendency can also be roughly seen in Fig. 7 of Zhang & 
MacFadyen (2009). 

Fig. 5 illustrates the effects of some parameters (n, 
-E'0,iso, ^j, and e) on the receding jet component in the af- 
terglow light curve. Fig. 5(a) shows that the circum-burst 
medium density (n) affects the peak time (tpoak) of receding 
jet dramatically. A larger number density usually leads to a 
smaller fpeak- The strength of the receding jet component is 
also obviously enhanced. It again hints that the receding jet 
component is most likely detectable in a dense environment. 
Similarly, the initial kinetic energy (-Eo,iso) also affects tpcak 
significantly, with larger -Bo,iso corresponding to a larger 
^peak (Fig- 5(b)). The effect of the initial jet opening angle 
(^j) on tpeak can also be clearly seen in Fig. 5(c). It should 
be further noted that the receding jet component is more 
marked when the opening angle is smaller. In Fig. 5(d), we 
can observe an obvious rcbrightcning when the radiation ef- 
ficiency (e) is large. However, in realistic case, e is unlikely 
to be so large. Actually, at such late stages, the external 
shock should be adiabatic, so that e should be nearly zero. 

In Fig. 5(d), we also plot the radio afterglow light curves 
for double-sided jets under some special physical assump- 
tions. The dash-dotted line is plotted by assuming that both 
the forward jet and the receding jet do not experience any 
lateral expansion. Since the deceleration of the jets is much 
slower in this case, we see that the receding jet component 



emerges much later and is also much less obvious as com- 
pared with our "standard" case. The dotted line is plotted 
by assuming a much smaller initial Lorentz factor (70 = 30), 
which may correspond to the so called failed GRBs (Huang 
et al. 2002). The receding jet component emerges slightly 
earlier as compared with the solid line, but its role becomes 
less significant correspondingly. 

Fig. 6 illustrates the effects of other four parameters (^c, 
^3, p, and ^obs) on the receding jet component. Generally 
speaking, a larger and/or can enhance the receding 
jet component markedly. On the other hand, although p has 
an important influence on the overall afterglow hght curve, 
its impact on the relative strength of the receding jet com- 
ponent is not significant. Again, note that in all the cases, 
the contribution from the receding jet only emerges as a 
plateau, but not as any obvious rebrightening. In Fig. 6(d), 
when the observing angle (f?obs) increases, the forward jet 
component becomes weaker, while the receding jet com- 
ponent becomes stronger. It is in good accord with our 
expectation (also see Granot & Loeb 2003). However, the 
contribution from the receding jet still generally plays a 
minor role in the total afterglow light curve. Additionally, 
for off-axis twin jets, the GRB from the forward jet is un- 
observable, so that even the afterglow from the forward jet 
itself (i.e., the orphan afterglow) is difficult to observe. Note 
that in Fig. 6(d), when ^obs = 7''/2 (i.e., the thick solid line), 
the contribution from the receding jet and the forward jet 
are actually equal. 

Equation (14) tells us that the peak time of the receding 
component should be relevant to the 3 parameters of n, 
£'0,iso, Oj] on the other hand, other parameters such as ^e, 

p do not affect the peak time. These tendency can be 
clearly seen in Figs. 5 and 6. 

In all the above calculations, we have assumed that the 
conditions and parameters of the twin jets are the same. 
However, this may not be the case for realistic GRBs. The 
circum-burst environment and the micro-physics parame- 
ters may actually be different for the twin jets, as that may 
happen in the two component jet structure (Huang et al. 
2004; Jin et al. 2007; Racusin et al. 2008). In Fig. 7, we 
have plotted the overall afterglow light curves by assuming 
different parameters for the forward jet and the receding 
jet. In each panel of Fig. 7, we first plot a common light 
curve (the solid line) by adopting the standard parameter 
set, but change to 0.01 and change to lO""*. We then 
increase the values of ^e, Ci' ^^'^ ^ ^he receding jet to 
see their effects on the afterglow light curve. It is encour- 
aging to see that the emission from the receding jet really 
can be greatly enhanced, so that it can manifest as an ob- 
vious rebrightening in the overall light curve. In Fig. 7(a), 
7(b) and 7(d), the peak flux of the rebrightening can be 
nearly 100 times larger than the "background" level in the 
best cases. It is imaginable that in the most favorable cases, 
when all ^c, Cb ^^'^ larger for the receding jet at the 
same time, the rebrightening will be even more remarkable. 
However, note that the contrary condition may also exist in 
realistic GRBs, i.e., these parameters may also be smaller 
for the receding jet. Then the emission from the receding 
jet will be completely unnoticeable. 

4. Conclusion and Discussion 

We have studied the emission of the receding jet numeri- 
cally. The effect of the EATS is included in our calculations. 
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Clearly, this effect plays an important role in the process. 
It is found that the contribution from the receding jet is 
generally quite weak. In most cases, it only manifests as a 
short plateau in the overall afterglow light curve, but not 
a marked rebrightening. The flux density of the plateau 
is usually much less than 100 /iJy in radio bands even at 
a small redshift of z = 0.1 . If we place the GRB at a 
more typical redshift of z = 1, then the flux density of the 
plateau will be less than 0.1 /zJy at 8.46 GHz. We noticed 
that the observed radio afterglow emission is generally on 
the level of 0.1 — 1 mJy at about the peak time. After 
several months, the radio afterglow usually decreases to a 
very low level, and is submersed by the emission from the 
host galaxy, whose strength can be 40 — 70 /xJy (Berger et 
al. 2001). Additionally, the error bar of radio observations 
is usually ~ 30 — 50 /xJy at very late stages (Frail et al. 
2003). Thus the contribution from the receding jet, i.e. the 
plateau, is actually very difficult to detect currently, espe- 
cially for those GRBs at z ~ 1. Our results are consistent 
with a recent observational report by van der Horst et al. 
(2008), who failed to detect any clear clues of the receding 
jet emission. 

However, as shown in our Fig. 7, if the micro-physics pa- 
rameters of the receding jet were different from the forward 
jet, or if the receding jet were in a much denser environ- 
ment, then it is still possible that the contribution from 
the receding jet can be greatly enhanced. For example, if 

and/or of the receding jet is much larger than that 
of the forward jet, then the receding jet can really manifest 
as an obvious rebrightening. 

Also, our Fig. 5(a) shows that a dense circum-burst en- 
vironment can suppress the emission of the forward jet, and 
enhance the contribution from the receding jet. If the GRB 
occurs in a very dense molecular cloud with n > 10^/cm^ 
(Dai & Lu 1999), the contribution from the receding jet 
may be much easier to detect. Additionally, if the GRB is 
very near to us at the same time, then the possibility of 
successfully detecting the receding jet is very high (see the 
thin lines in Fig. 3(a)). 

In short, we believe that the effort of trying to search 
for the afterglow contribution from the receding jet is still 
meaningful. If observed, it would provide useful clues to 
study the circum-burst environment and the micro-physics 
of external shocks. We suggest that nearby GRBs (with 
redshift z < 0.1) should be good candidates for such stud- 
ies. 
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Fig. 1. The evolution of the Lorentz factors of the twin jets. The sohd hne corresponds to the receding jet and the dashed 
hne is plotted for the forward jet. The twin jets are in "standard" condition as defined in Sect. 3. The observers' time 
has been corrected for the cosmological effect {z = 0.1). 
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Fig. 2. Schematic illustration of the EATSes at three moments, ti « 114 d, t2 « 286 d and ts ~ 572 d. In this calculation, 
we have used the "standard" parameter set as defined in Sect. 3. "O" is the position of the central engine, and the observer 
is on the far right side with Y=0. The dotted lines indicate the jet boundary. For the receding jet, the EATSes are plotted 
in thick solid lines, while for the forward jet the surfaces are plotted in thin solid lines. Note that on the forward jet 
branches, the bulk Lorentz factors of the material at the peak of the EATSes arc 1.17, 1.07, and 1.03 for ti, t2, and t^, 
respectively. On the receding jet branches, the bulk Lorentz factors of the material at the peak of the EATSes are 56.07, 
11.79, and 3.95 for ti, t2, and ts, respectively. 
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Fig. 3. 8.46 GHz radio afterglow (a) and R-band optical afterglow (b) from the forward jet and the receding jet. The thick 
lines are plotted for a "standard" double-sided jet as defined in Sect. 3. The thin lines are plotted for the double-sided 
jet with only one parameter altered as compared with the "standard" condition, i.e. n = 1000/cm^. In each group, the 
dotted line reflects the emission from the forward jet, the dashed line reflects the contribution from the receding jet, and 
the solid line is the total light curve. 




Fig. 4. Multiwavelength afterglow light curves of a double-sided jet. Radio afterglows are illustrated in panel (a), and 
optical/IR afterglows are plotted in panel (b). In this calculation, we have used the "standard" parameter set as defined 
in Sect. 3. 
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Fig. 5. The effects of various parameters (n, -Eo.iso, Sj, and e) on the 8.46 GHz radio afterglow hght curves of double- 
sided jets. In each panel, "(s)" corresponds to the "standard" condition as defined in Sect. 3, and other lines are drawn 
with only one certain parameter altered or one condition changed. In panel (d), the dash-dotted line is plotted for a 
double-sided jet without lateral expansion; and the dotted line is plotted for a double-sided jet with a low initial Lorentz 
factor (70 = 30), which may correspond to the so called failed GRBs. 
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Fig. 6. The effects of various parameters (^ej Cb! P> ^'^^ ^obs) on the 8.46 GHz radio afterglow hght curves of double-sided 
jets. In each panel, "(s)" corresponds to the "standard" condition as defined in Sect. 3, and other lines are drawn with 
only one certain parameter altered. 
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Fig. 7. 8.46 GHz radio afterglow light curves of double-sided jets. In this figure, we assume that the parameters of the 
receding jet can be different from those of the forward jet. In each panel, the solid line is plotted under the "standard" 
condition, i.e., the parameters arc completely the same for the twin jets (but note that wc have evaluated as 0.01 and 
as 10~^ here). For other light curves, one or two parameters are changed for the receding jet, to see its effect on the 
afterglows. 



